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Abstract 

We explore the feasibility of detecting heavy neutrinos by the existing facilities of 
neutrino experiments. A heavy neutrino in the mass range 1 MeV < Mn Hz 500 MeV 
is produced by pion or kaon decay, and decays to charged particles which leave 
signals in neutrino detectors. Taking the T2K experiment as a typical example, we 
estimate the heavy neutrino flux produced in the neutrino beam line. Due to massive 
nature of the heavy neutrino, the spectrum of the heavy neutrino is significantly 
different from that of the ordinary neutrinos. While the ordinary neutrinos are 
emitted to various directions in the laboratory frame due to their tiny masses, the 
heavy neutrinos tend to be emitted to the forward directions and frequently hit 
the detector. The sensitivity for the mixing parameters is studied by evaluating 
the number of signal events in the near detector ND280. The sensitivity of T2K 
at 10 21 POT is found to be comparable to that of the previous experiment PS191, 
which has placed the most stringent bounds on the mixing parameters of the heavy 
neutrinos for 140 MeV < M N < 500 MeV. 



1 Introduction 



In the last few decades, neutrino oscillation experiments have conclusively shown that 
neutrinos are massive pQ. The minimal version of the Standard Model is thus to be ex- 
tended to accommodate the neutrino masses. In many possible extension of the Standard 
Model, neutral heavy leptons are often predicted. In the seesaw mechanism [2] for ex- 
ample, the right-handed neutrinos are introduced and they weakly mix with the ordinary 
neutrinos after the electroweak symmetry breaking. 

For the mass of the heavy neutrinos, a wide range of possibilities has been discussed 
in literature. In the canonical picture of the seesaw mechanism, heavy neutrino masses 
are supposed to be around the grand unification scale. These super heavy neutrinos 
can account for the baryon asymmetry of the universe by the leptogenesis scenario [3J. 
Another possibility to account for the baryon asymmetry has been suggested in jU |S] 
and further studied in [6]. In this scenario, two quasi-degenerate heavy neutrinos of 
(9(100) MeV — O(10) GeV play a crucial role in the early universe. Heavy neutrinos in the 
mass range ~ 0.2 GeV could enhance the energy transport from the core to the stalled 
shock and favor the supernova explosion [7] . Heavy neutrinos with a few keV mass have 
also attracted much interests as a viable dark matter candidate [S] and an agent of the 
pulsar velocities [9]. Remarkably, the dark matter and the baryon asymmetry due to keV 
and GeV heavy neutrinos can originate in a simple framework so called z/MSM [TO], [5] , 
which is an extension of the Standard Model with just three generations of the right- 
handed neutrinos. Besides the super heavy range much larger than TeV, such heavy 
neutrinos can be tested in existing and forthcoming experiments due to lower threshold 
energies of production (for example, see Ref. [HI [12] and references therein). 

In this paper, we focus on the heavy neutrinos produced by kaon decays. The previous 
neutrino experiments, including peak searches in the meson decays [T3] Q31 [151 EEl 02] an d 
the decay search with accelerators [TBI [T9] , have placed stringent bounds on the mixing 
parameters in this mass range. In particular, PS191 [19] has placed the strongest bounds 
for the mixing parameters in the mass range of 140 MeV < < 500 MeV. Since the 
PS191 experiment in 1984, however, no further experiments of this type of decay search 
have performed and the bounds have not been updated for about 30 years. On the 
other hand, great progress has been made in neutrino oscillation experiments over the 
same period. It is interesting to note that typical long or short baseline experiments are 
equipped with the (near) detectors placed at 0(100) meters away from the beam targets, 
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T2K [23J 


MINOS [24J 


MiniBooNE [25 


SciBooNE [26J 


POT 


0.86 x 10 19 


10 21 


10 21 


10 21 


10 21 


(Distance) -2 


(128 m)" 2 


(280 m)" 2 


(lkrn)" 2 


(541m)- 2 


(100 m)- 2 


Volume 


216 m 3 


88 m 3 


303 m 3 


524 m 3 


15.3 m 3 


Events 


1 


9.9 


2.7 


15.8 


13.5 



Table 1: A comparison between PS191 and recent accelerator experiments. The item 
"Distance" means the distance between the beam target and the detector for each ex- 
periment. The item "Events" shows POT x (Distance) _2 x Volume in units of PS191. The 
POTs for the oscillation experiments are assumed to achieve 10 21 . 



which detectors are capable of measuring the charged-particle tracks produced by the 
heavy neutrino decays. A natural question is then whether the existing accelerator-based 
neutrino experiments are capable of discovering the heavy neutrinos and how sensitive 
such experiments arqj. We believe that this is a timely question to ask. In fact, the 
exposure of PS191 is about 200 m 3 x 10 20 POT while the recent accelerator-based neutrino 
experiments are expected to achieve 10 21 POT with the near detectors which are typically 
no smaller than 200 m 3 by factor of 10. Table [1] shows a comparison between PS191 and 
several examples of recent neutrino experiments. Among several options of accelerator 
experiments, in this paper we focus on the T2K experiment as a typical example. 

The study follows two main steps; the flux estimation and the event number calculation 
for various signal decays. First, in the fulx estimation, we use a semi-analytical method 
with a help of the active neutrino flux <j) l/ simulated by the T2K collaboration [27] [28| 129] . 
making a reasonable simplification for the geometry of the decay tunnel and the detector. 
More precise fluxes of the heavy neutrino might be obtained by Monte Carlo methods. In 
a Monte Carlo calculation it is possible to take into account the details of the geometry 
and the spectrum of the parent mesons. The analytical technique is nevertheless useful, 
since it allows one to understand the essential physics which determines the behavior of 
the heavy neutrino spectrum, in particular its mass dependence. We emphasize that the 
phase-space effect in kaon decay is important and the heavy neutrino flux (j)^ can be 
significantly deviated from the naive expectation <Pn — |©| 2 0^, where O is the active- 
heavy mixing parameter of interest. As the heavy neutrino mass Mjv approaches to the 
production threshold, the heavy neutrinos tend to be distributed to lower energies with a 
narrow spread, so that the proportionality 4>n — |©| 2 0f is broken. This phase-space effect 



*Far detectors with large volume are also useful to detect the heavy neutrinos produced in the atmo- 
sphere [201 HI] . 
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leads to larger event numbers than the naive expectation and has a significant impact for 
the estimation of the sensitivity. 

Second, in the event number calculation, we take into account various decay modes 
of the heavy neutrino N. The two-body modes N — y e T 7r ± and N — y / u =F 7r ± are the most 
promising channels due to their large branching ratios. For these modes, the invariant 
mass distribution for the lepton and the pion momenta has a peak at the heavy neutrino 
mass. The three-body modes N — y e~e + v, N —y fi T e ± u and iV —y fi~fi + u are also 
interesting. While the first mode may suffer higher background by tt° decay, the latter 
two modes have smaller backgrounds and also serve as promising channels for discovering 
the heavy neutrino. By assuming non-observation of these exotic events in the T2K near 
detector at 10 21 POT, one finds the upper-bounds for the mixing parameters comparable 
to PS 191. This means that T2K may have a good chance to discover the heavy neutrino 
in near future. 

The layout of this paper is as follows. In Section 2, we introduce the heavy neutrino 
and briefly review its essential properties. In Section 3, the flux of the heavy neutrino 
at the near detector is discussed. In Section 4, the decays of the heavy neutrino at the 
detector and their event numbers are studied. Section 5 is devoted to conclusions. 

2 Properties of the heavy neutrino 

We consider a heavy (sterile) neutrino N with the mass range 1 MeV < M N < 500 MeV. 
The flavor eigenstates of the left-handed neutrinos v a (a = e, //, r) are given by the linear 
combination of the mass eigenstates as 

v a = U ai Ui + Q a N, (1) 

where U ai is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, z/j (i = 1,2,3) 
are the mass eigenstates of the active neutrinos. The parameter Q a is the mixing be- 
tween (light) active and sterile neutrinos, through which N interacts with the weak gauge 
bosons. The extension to the multi-generation case is trivially done by replacing O a iV 
with J2i ®aiNi. In this paper, we assume iV is Dirac particle unless otherwise stated. 
This is simply because we would like to make a comparison to PS191 in which the same 
assumption is made. If N is Majorana particle, the decay width is doubled since N decays 
to charge-conjugated states also. 



3 



With the mass range of 140 MeV < Mn < 500 MeV, the heavy neutrino N is produced 
by K decay. The main production modes are 

K + fi + N, K + e + N (2) 

for M N < 388 MeV and M N < 493 MeV respectively. The decay width K + fi+N 
(K + — > e + N) is proportional to | ©^1 2 (|0 e | 2 ), and G r is irrelevant for the production^. 
We describe further details of the decay modes (T5]) in Section 3 and Appendix A. Since the 
magnetic horn focuses positively-charged mesons, the contributions from K~ are small [28] 
and in what follows we neglect the K~ contributions. 

The heavy neutrinos N produced by the meson decays escape from the decay volume 
and part of them are injected into the near detector ND280 in the T2K experiment and 
decay to leave signals. Depending on the mass, N decays to lighter particles in various 
decay modes; 

N ->■ 7z/, N -> 3u, N e~e + z/, N ->■ fi T e ± u, 

With the ND280 detector capable of identifying e^, ^ and 7r ± , some of the above channels 
can be detected. 

Due to large branching ratios, the two body modes N — > un , N — > e~n + and N — > 
[i~7i + would be the most frequent events. However iV — > utt° does not seems to be as 
promising as the other two modes since 7r° are copiously produced by the ordinary neutrino 
interactions and may suffer from large backgrounds. On the other hand, iV — > e~7t + and 
N — > fi~7i + leave two charged-particle tracks with monochromatic energies in the rest 
frame of N. The invariant mass distribution of the two charged-particle momenta sharply 
peaks at the heavy neutrino mass, which seems to be easier to see than a slight excess 
of 7T° events over huge background. We thus proceed without any further analysis of 
N — > wit , but study N — > e~7i + and N — > /i~7r + in more detail in Section 4. The decay 
rates of N — > e~7r + and iV — > fi~ir + are proportional to |G e | 2 and |S/J 2 respectively. The 
radiative decay N — > jv is induced at one loop and negligible in this work. 

In spite of small branching ratios, the three body modes iV — > e~e + u, N — > / u =F e ± z/ and 
N — > are also interesting to study. Especially N — > /i =F e ± z/ and iV — > /i + z/ have 

smaller background than iV — > e~e + v and a few these events may lead to the discovery of 



^The heavier mesons such as D can decay producing tau so that Q T is involved in the production. In 
this work we neglect the contribution from such heavier mesons. 
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Figure 1: Summary of the production and the detection processes of the heavy neutrino 
JV. The decay mode iV — > vrr is abbreviated here because it is not as promising detection 
channel as the other ones. 

the heavy neutrino. Notice that iV — > e~e + i> and iV — > pTp^v are conducted not only by 
the charged-current (CC) but also by the neutral current (NC), so that these decay rates 
depend on all types of the mixing parameter Q e ,fi,r- On the other hand N — > / u T e ± z/ are 
mediated only by the CC interactions and the decay rate depends on ei/i . 

Fig. [1] summarizes the production and the decay processes of N to be studied in this 
work. We do not present the formulas for the decay width here. A complete list of the 
decay modes and the decay width are found in Ref. [TT] . 

3 Heavy neutrino flux at the near detector 

In the T2K experiment, pions and kaons are produced by the interaction of 31 GeV protons 
with the graphite target. The produced mesons are focused by the magnetic horns and 
enter the decay volume of 96 m long filled with helium gas. The parent mesons decay in 
flight inside the decay volume. The off-axis detector ND280 is located 280 m from the 
target station. The off-axis angle to ND280 from the target position is 2.04°. The layout 
of the secondary beam line and the near detector is sketched in Fig. [2j Details of the 
experiment setup are found in Ref. [29] . 
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Beam dump 



Figure 2: Schematic of the secondary beam line and the near detector ND280. 

Calculation of the neutrino flux at the near detector is a complicated task. The T2K 
collaboration has simulated the fluxes of the active neutrinos by the Monte-Carlo method. 
In this work, we do not follow their approach, but estimate the heavy neutrino flux by 
a semi-analytical method similar to Ref. [22]. We try to reconstruct a reasonable flux of 
parent particles from the active neutrino flux of Ref. [28], and then evaluate the heavy 
neutrino flux from the reconstructed parent flux. In this paper, we focus on K + meson 
as the parent of the heavy neutrinos since K + decay covers wider range of the heavy 
neutrino mass than 7r + decay. The following discussion is however easily extended to the 
7r + case. 

3.1 Modeling the parent flux 

For ND280, the neutrino source is a line-like object rather than a point. Let 4>(pk, 
denotes the K + spectrum along the decay volume, where px is the magnitude of the K + 
momentum and I is the flight length of K + . We set / = to be the upstream end of 
the decay volume. In the decay volume filled with helium gas, the decay length of K + is 
much shorter than the interaction length. One thus finds 

i_ 

<Pk(pk,1) = §k{pkY a k , (3) 

where 4>k{pk) is the spectrum at / = and Ak is the total decay length Ak = ^^{pk/^k) m 
with the kaon mass ttik = 493 MeV. 

If the parent spectrum 4>k{pk) is known, one can calculate the daughter flux from 
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K + decay. The source term of is given by 

(I>k(pk,i) i ^ 3 r 



S v (E v ,8,<f>,l) = / dp K 
Jo 



/3(l/r) 7 TdE u d cos 



^, /f(p/ , (^)_i5_, (4 ) 



where E v is energy, 9 and <fi are the polar and the azimuth angles of the emitted 
relative to the K + momentum directions, and V is the K + — > /U + z/ M decay width. Giving 
the source term, the z/ M flux <p Ufi (E u ) at ND280 is obtained by integrating the source term 
along I and the angles 9 and <fi covered by ND280. It reads 

<j> Vlt (E v ) = Jdl J dcos9 d<P^S v (E v ,9,<p,l)P(9,<p), (5) 

where If = 96 m, A is the effective area of ND280, P(9,4>) is the "projection" function 
which turns unity if enters ND280 and otherwise zero. 

The parent kaons carry two more degrees of freedom which are not explicitly mentioned 
above; the polar angle relative to the beam axis (the central axis of the decay volume) 
and the radial coordinate from the beam axis, both of them are defined at I = 0. The 
location of ND280 relative to each K + momentum depends on these degrees of freedom 
and the function P(9, 0) is in fact highly complicated. However the situation is greatly 
simplified if kaon momenta are assumed to be parallel to the beam axis. Information 
about the polar angle and the radial coordinate is then to be represented by an effective 
off- axis angle 9q, which is an "virtual" off-axis angle to ND280 from the upstream end of 
the decay volume. This angle #o represents the average off-axis angle to ND280 from each 
K + momentum, which angle varies kaon by kaon carrying deferent polar angles and radial 
coordinates. Furthermore we neglect 9 dependence of the effective area A. Adopting these 
simple modeling of geometry of the decay volume and the detector, we use 

<j> Vll {E v ) = ^ Jdl J dcos9S u (E u ,9,l)P\9,9 ), (6) 

as a formula to relate the daughter z/ M flux with the parent K + spectrum. Here A0 
is cf) interval defined by the detector width, P'(9, 9q) denotes the projection function 
determined by the detector hight and K + decay point /. Since the explicite expression of 
Eq. is rather lengthy, we put details on Eq. fl6]) in Appendix |A] . 

Having Eq. ([6]), our strategy is fitting <p v ^{E v ) calculated in Ref. [28] by adjusting 
the parameters included in the right-handed side of Eq. ([6]), with a kaon spectrum which 
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E v (GeV) 



Figure 3: A comparison between Eq. fl6]) and the flux simulated in Ref. [28J. Dots shows 
the result of the simulation in Ref. [28] and solid curve is Eq. ([6]) with the parameters 
# = 1-48°, a = 4.8 x 10 19 mb" 1 and p = 2.1 GeV. 

is physically well-motivated. The K + spectrum da/dpx in the proton collision with a 
graphite target is measured by NA61/SHINE Collaboration [30], which is customized to 
improve the flux calculation in T2K. It seems reasonable to expect that the shape of the 
kaon spectrum is not far from this measured spectrum. However, the effects such as the 
secondary protons [2B] and the magnetic horns must deform da/dpx to some extent. To 
take into account this deformation, we allow a shift of the peak of da/dpx- In summary, 
in order to model the kaon spectrum 4>k{pk), we introduce two free parameter ao and po, 
the overall scale factor and the shift of the peak, respectively. See Appendix |A] for more 
details. 

With the above ansatz for the kaon spectrum, we have three parameters in the right- 
hand side of Eq. ([6]); the effective off-axis angle 6q, the scale factor ao and the shift of 
the peak pq. Fig. [3] shows the fit by Eq. ([6]). The dots show the result of the simulation 
in Ref. [28], while the solid curve shows Eq. ([6]) with the parameter set of 6$ = 1.48°, 
a = 4.8 x 10 19 mb _1 and po = 2.1 GeV. It is seen that Eq. (jSJ) is well tracing the global 
behavior of the simulated result of Ref. [28] . 

Notice that the effective off-axis angle is taken as 6$ = 1.48°, which is smaller than 
the actual one 2.04°. If this angle is taken as #o = 2.04°, the flux starts to fall off at 
E v ~ 6 GeV and does not fit the data points above 6 GeV. This is simply due to the 
two-body kinematics. In K + — > decay, the maximum energy of is 

2 2 

mi — mt 

v ~2m^sin# K} 
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Figure 4: Fluxes of the heavy neutrino <Pn{En) from the K + — > fi + N mode for several 
sample values of M^- Black dots shows jG^I 2 ^ • The masses are taken as = 
350 MeV (purple), 300 MeV (blue), 200 MeV (cyan) ^100 MeV (green), 50 (orange) MeV. 

for 9 < 90°. The larger the off-axis, the smaller the maximum neutrino energy. The 
abundance of above 6 GeV means that certain fraction of kaons must have momentum 
directing to the near detector so that with smaller 9 can contribute to the flux. This 
is of course expected since the magnetic horns do not make the kaon momenta perfectly 
parallel to the beam axis. The smaller angle 9q = 1.48° effectively takes this effect into 
account. 

3.2 The heavy neutrino flux 

With the kaon spectrum <Pk(pk) discussed above, let us estimate the heavy neutrino flux 
4>n(E]s[). The calculation goes as before. The source term is given by 

where Ti and T 2 are the decay width for K + — > fi + N and K + — > e + N respectively. 
Provided that the decay length of the heavy neutrino is much larger than the distance 
between the K + decay point and the detector, the heavy neutrino flux <Pn(En) is given 
by just replacing S u with in Eq. ([6]). We will see in Section 0] that this is the case for 
the energy and the parameters of the current interests. Details on the differential decay 
rates dTi are shown in Appendix lAl 

Fig. H] and [5] present the heavy neutrino fluxes (()n{En) for several fixed values of M^. 
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Figure 5: Same as Fig. H]but for the K + — > e + N mode. 

Fig. H and [5] show the case where (|0 e |, |9 M |) = (0, 10" 1 ) and (|O e |, |0 M |) = (KT^O) re- 
spectively. The masses are sampled as M N = 350 MeV (purple), 300 MeV (blue), 200 MeV 
(cyan), 100 MeV (green), 50 (orange) MeV. In the both cases of K + — > fi + N and K + — > 
e + N, the spectral shapes are similar to 0^ for < 100 MeV, whereas they are sig- 
nificantly deviated from (j) u for Mn > 100 MeV. The most remarkable feature is the 
enhancement of the flux at lower energies for larger M^. In Fig. HJ for example, the 
heavy neutrinos of Mn = 350 MeV gather around 2 — 3 GeV and the peak intensity 
reaches 10 11 cm _2 GeV _1 . This is nearly two orders of magnitude larger than the naive 
expectation |6 M | 2 0^ M . 

There are two reasons for the enhancement. One is the spin conservation in the K + 
rest frame. The matrix elements of the decay processes K + — > fi + N and K + — > e + N scale 
with M^r when Mn is much larger than the mass of the charged lepton in each final state. 
This accounts for the larger fluxes than \Q e ^\ 2 (t) Vll for > m^. This also accounts for 
the smaller flux than |0 e | 2 0^ M for M N = 50 MeV in Fig. |5j 

The another reason is slower motions of the heavy neutrinos at the rest frame of the 
parent particle. The smaller the daughter velocities, the easier to boost them into the 
forward directions. In the current setup, the detector may well be regarded as the object 
placed in the forward direction of the K + momenta. In the case of K + — > /i + z/^, neutrinos 
can be emitted not only to the forward directions but also to the backward directions since 
the neutrino masses are very small and the neutrino velocities at the rest frame are larger 
than the parent typical velocity at the laboratory frame. On the other hand, in the decay 



10 



K + fi + N, the heavy neutrinos tend to be emitted to the forward directions. In the 
rest frame of K + , the gamma factor of N at the K + rest frame is given by 

= <Z<±^L . (9) 

This goes to unity as approaches to the threshold value = mx — ra^ = 388 MeV. 
On the other hand, most of the kaons carry the momentum around 1 — 4 GeV, so that the 
gamma factor of K + is typically given by jk = 2 — 8. Then Eq. (jUJ) tells us that jk > 7jv 
for Mn > 120 MeV. Hence for > 120 MeV, kaon's velocities overcome N's velocities 
and N are focused into the forward directions. This agrees with the flux behavior seen in 
Fig. land EJ 

4 Event rates and expected sensitivity 

As we have seen in Section [2J a fraction of the heavy neutrinos passing through the 
detector decay inside the detector volume and leave signals via various decay modes. In 
this section, we calculate the number of signal events in ND280 and estimate the potential 
sensitivity of ND280. We argue that the sensitivity of ND280 is comparable to that of 
the PS 191 experiment. 

4.1 Number of the signal events 

The total number of events is given by the difference between the number of the heavy 
neutrinos at the up and down stream end of the detector. For a particular decay channel, 
the number of events is given by 

i>oo ^ j-X 1 

Events = A I dE N — I dx 4>n{E N) x), (10) 

J M N * Jx 

where A is the (partial) decay length for the signal decay mode of interest, A is the cross- 
sectional area of the detector, x is the flight distance of the heavy neutrinos and (xq,xi) 
means the detector segment. The number of the heavy neutrinos decreases by the decays. 
With the total decay length Ajv, the x dependence of 4>n{En,x) is determined as 

(f) N (E N ,x) = <p N {E N )e~^, (11) 

where (Pn(En) is the heavy neutrino flux discussed in Section I3T21 
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M N (GeV) M N (GeV) 

Figure 6: Parameter region that satisfy A^r < 300 m for 7^ = (The filled region 
labeled by A^v < 300m). The filled region labeled by r > 0.1(s) shows the region where 
heavy neutrino decay may spoil BBN. 

Eq. ( TlOl is further simplified if the total decay length is much larger than the flight 
distance of the heavy neutrino. Provided that A a? ^> xq,x\ — Xo, the number of events 
reads 

r°° v 

Events ~ / dE N (j) N (E N ) — , (12) 

JM N * 

where V = A(x\ — Xq) is the detector volume. In the T2K setup, xq ~ 300 m and 
ii-io~5m. Thus the condition Ajv 3> X\ — xq holds if A^ ^> x 300m. 

It turns out that the condition 3> 300 m holds in the most parameter and energy 
region of interest. Fig. [6] highlights a "strong coupling" regime where Ajv < 300 m. Here 
the region is showing an example with 7^ = \/2 (pn = Mat)- For 7at = a/L5 (pn = 
Mjv/2), the boundary is pushed down by factor of two. For the energies 7^ > y/2, the 
total decay is effective only for strong coupling regime | © e ,^t 1 2 ^ 10~ 4 , which is already 
ruled out by many experiments. Eq. ffl2|) is well applicable for the following discussions 
which treat the region far below the bound shown in Fig. |6l 

In Fig. |6l we put in passing the region where the lifetime of the heavy neutrino becomes 
long enough so that late time decay of the heavy neutrinos may spoil the success of 
Big Bang Nucleosynthesis (BBN). Ref. [311 [32] has studied such a bound for lOMeV < 
Mm < 140 MeV in detail. For Mn > 140 MeV, however, there is no consensus about the 
constraint from BBN. Here we simply present the region for r > 0.1 s [32] where the heavy 
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Figure 7: Expected sensitivity of T2K and upper limits by PS191. The blue solid curves 
show the contour for 2.44 events at T2K. The red curves show 90% CL upper bound by 
PS191 ng. 

neutrinos are not cleared away before the onset of BBN. 
4.2 Comparison between T2K and PS191 

For each mass eigenstate of the heavy neutrinos, four additional parameters are intro- 
duced into the Standard Model; and B e>jU>T . Experiments impose some constrains on 
the four dimensional space (Mjv, |@e|, I ©/J, |©r|)- However, the analysis of experimental 
constraints on the full four- dimensional space is a complicated task. PS 191 has made the 
following assumptions/simplifications in their analysis. 

• Heavy neutrinos are Dirac particles. 

• The NC contributions to the three-body decays of iV are neglected. 

• Either K + — > fi~N or K + — > e~N is dominant in the production. 

In the following, we first make the same simplifications, just aiming for a comparison 
between T2K and PS191. In Section 14.31 we give a comment on the case where the 
second simplification is relaxed so that the three-body decay depends on T . 

4.2.1 Two body channels 

Let us first focus on the two body decays. Fig. [7] shows the expected sensitivity for 
the chains K + -»■ e + N -»■ e + (e~7r + ) (left), K+ e + N e + (/i~7r + ) (middle) and 



13 



K + ^e + N^e + (e-e + v e ) K*-*ifN-+tf{*-0*vJ ^ + ->/* + N->/* + (/re + v e ) 




0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 

M N (GeV) M N (GeV) M N (GeV) 

Figure 8: Same as Fig. [7] but for the three-body decay channels. 



K + — > [i + N — > /i + (//~7r + ) (right), respectively. In the figures, the red curves show 90% CL 
upper bound by PS 191 [19]. The blue solid curves show the contour for 2.44 events, which 
corresponds to 90% CL limit when the measured signal and the expected background are 
null [33J. Here the fiducial volume of the detector is taken as V = 3.5 x 3.5 x 5.0 = 
61.25m 3 [31]. 

The interactions between the active neutrinos and the nuclei in the detector provide 
backgrounds for the decay signals N — > h~tt + and iV — > e _ 7r + . For instance, the reactions 

+ n — > A*~ + 7r+ + n (CC — nn + ) 

Vy, + 16 fi- + 7T+ + 16 (CC - coherent tt+) 

may become background for iV — > fi~n + . It is expected that these events account for 4% 
of the whole neutrino events in ND280, resulting 7300 events/10 21 POT/ton [31]. However 
the background can be reduced by taking the invariant mass of fi~ and 7r + momenta in 
the final state. For the heavy neutrino signal, the event distribution sharply peaks at 
the heavy neutrino mass while the events provide continuous background. A serious 
sensitivity should be estimated together with the invariant mass distribution for the 
reactions, the energy resolutions, all sort of uncertainties, etc. Such a thorough analysis 
is interesting but beyond the scope of this work. 

We can further reduce the background by selecting the events taking place in the TPC 
volume which is filled by argon gas. Due to low density of the gas region, the events are 
significantly reduced while keeping the signal rates unchanged. Out of the full volume of 
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61.25 m 3 , 9.0m 3 is filled with the argon gas [34] and available for this purpose. According 
to Ref. [35], the total neutrino events taking place in the argon gas are about 2000 at 
10 21 POT. Since 4% of them become the background, the number of the background event 
is expected to be around 80. These 80 events will be further reduced in the bin around 
the heavy neutrino mass. Keeping this in mind, we plot the contour for 2.44 events with 
V = 9.0 m 3 by the dashed curves. 

For iV — > e~iT + , the background processes are produced by the CC interactions of v e . 
However, the v e flux is about two orders of magnitude smaller than that of around 
the peak energy ~ 600 MeV [271 ESI EH]- By selecting the events in the gas volume, 
the background rate is expected to be less than one for 10 21 POT. The decay channel 
iV — > e~7r + is more promising than N — > ^~tt + in view of signal/background ratio if 

|e e |~|e M |. 

4.2.2 Three body channels 

Fig. [8] presents the same plots as Fig. [7] but for the three-body channels; K + — > e + N — > 
e+(e-e + z/ e ) (left), K + -> /j,+N ^ + (e~e + v e ) (middle), K+ fi+N fi + (fi~ e+ v e ) 
(right). 

A major obstacle to successful identification of iV — > e~e + u would be tt° that is 
copiously produced by neutrino interactions. The two photons from tt° decay develop 
to electromagnetic cascades in the detector material and may mimic the signals. The 
subdominant decay mode 7r° — > e~e + 7(1.17%) may also contribute to the background 
when one of the final-state particle is undetected. The invariant mass distribution of 
the electron pair is useful since it moderately peaks at one-half of the heavy neutrino 
mass |21j. The analysis needs anyway precise understanding of the background, and the 
detection via iV — > e~e + v seems less promising than the two-body modes. 

As for iV — > /i~e + z/, the charmed- meson production by the neutrino CC interaction [36] 
and successive semi-leptonic decay may become the background. According to Ref. [36] . 
the cross section of the charm production is about 1% (4%) of the total CC cross section 
for E v = 5GeV(15GeV). Due to the off-axis technic, however, the contributions from 
such high-energy neutrinos are suppressed in the T2K setup. By selecting the events in 
argon gas, the neutrino reduction rates can be further reduced, and iV — > \T e + v may 
become more or less background free. 

Although PS191 has not studied the signal process iV — > fi~ open for M N > 
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t>0.1(s), 10,1=10^1=0 
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Figure 9: Allowed region of Mjv-|9 r | 2 plane. The upper- filled regions are excluded by 
CHARM [37] and DELPHI [HE] at 90% and 95% CL, respectively. The lower-filled region 
is the regime where the lifetime of the heavy neutrino is longer than 0.1 s. 

211 MeV, it should be emphasized that searching for this di-muon signal is also a promising 
method. The main background for iV — > {i~fi + v may be the charmed-meson production 
by the neutrino CC interaction [36] and successive semi-leptonic decay as in the case of 
N —7- [i~e + v. This rate is however expected to be small for the neutrino energy in T2K. 

4.3 Implications of G T ^ 

So far we have focused on the comparison between T2K and PS191 and T is accordingly 
neglected. In this subsection, we comment on several implications of the G T ^ case. 

Since G r is not involved in the main production processes such as pion and kaon 
decays, the experimental constraints of |0 T | 2 are much wearker than that of |0 e | 2 and 
\Q^\ 2 . Fig. IH] shows the allowed region of Mtv-|0 t | 2 plane. The upper-filled regions are 
excluded by CHARM [37] and DELPHI [38J at 90% and 95% CL, respectively. The lower- 
filled region is the regime where the lifetime of the heavy neutrino is longer than 0.1s. 
By comparing Fig. [9] and Fig. [HI it is seen that, unlike |O eiA t| 2 , |@r| 2 can be large enough 
so that the BBN constraint is avoided without contradicting with the upper bound from 
the direct experiments. 

While 6 T is not involved in the production processes, the detection processes in general 




16 



c±_ 10" 5 




© 10- 9 — 

10 -io| . 1 . 1 . 1 

0.20 0.25 0.30 0.35 0.40 0.45 0.50 

M N (GeV) 



Figure 10: Expected sensitivity for N — > fi |© M | = is assumed. 

depend on T since N — > e~e + v and iV — >■ \T p + v are induced not only by the CC but also 
by the NC interactions (20] . This means the experimental setup discussed in this paper has 
sensitivities to certain combinations of |© e ,^| 2 and |© T | 2 - To demonstrate this, let us focus 
on a simple case where |© e | 2 , |© T | 2 ^ |©/t| 2 an d Mat > 211 MeV. In this case, the heavy 
neutrinos are produced by K + — » e + N and can be detected by N — > pT ji + u. Since l©^! 2 
is small, the decay process N — > fi~fi + i> is conducted only by the NC interactions and the 
decay width becomes proportional to |© e | 2 + |@r| 2 - Therefore by analyzing the dimuon 
signals in the detector, one is able to constraint |6 e |^|6 e | 2 + |6 r | 2 for M N > 211 MeV 
(or discover the heavy neutrino). 

Fig. [10] shows the expected sensitivity for iV — > As in Fig. [7] and Fig. [8] 

the solid (dashed) curve is the contour for 2.44 events with V = 61.25 (9.0) m 3 . From 
Fig. [10] and Fig. [9] one can see that there exists the parameter regime where the dimuon 
signal is sizable while the success of BBN is unspoiled. For example, it is seen that 
the combination |6 e | = 10" 4 ' 5 and |6 r |- 2 - 5 is BBN safe but O(10) - C(10 2 ) events of 
N — > fi~ [i + v are expected for 300 MeV < M N < 400 MeV. Interestingly, this case also 
predicts 0(1) — 0(10) events of N — > e~ir + (see the left pannel of Fig. [7]), so that the 
heavy neutrino model can be tested in multi-dimensional way. 
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5 Conclusions 



In this paper, we have focused on the heavy (sterile) neutrinos produced by kaon decays 
and explored the feasibility of their detection at the existing facilities of the accelerator- 
based neutrino experiment. Taking the T2K experiment as a typical example, we have 
estimated the heavy neutrino fluxes produced in the beam line and calculated the event 
rates of their decay taking place inside the detector. 

Due to massive nature of the heavy neutrino, the spectrum of the heavy neutrino is 
significantly different from that of the ordinary neutrinos. The ordinary neutrinos are 
emitted to various directions in the laboratory frame due to their tiny masses. On the 
other hand the heavy neutrinos carrying a large mass tend to be emitted to the forward 
directions and frequently hit the detector. This is a unique advantage of the experiments 
in which the parent mesons decay in flight with sufficient gamma factors. 

Among various decay modes, N — y e~n + open for > 140 MeV is one of the most 
promising channels for detection because of its larger rate and lower background. The 
backgrounds from the active neutrino reactions can be reduced by selecting the events 
occurring in the regions filled with no material. In the T2K near dector, the TPC volume 
9 m 3 filled with argon gas plays this role. The expected sensitivity for this mode is 
comparable to that of PS 191, which has placed the most stringent bound on the heavy 
neutrino mixing. 

The three body modes N — y e~e + u, N — y fi~e + v, and iV — y {i~fi + v are also interesting 
signals to search for. In particular, iV — y e~e + v and iV — y [i + v are conducted not only 
by the charged current but also the neutral current, so that the tau flavor mixing G T is 
involved in the detection probabilities. Since |B r | is less constrained than |0 e | and |0 M |, 
the above two modes are not necessarily suppressed when |G e | and |G M | are small such 
that the two-body modes N — y e~7r + and iV — y fi~7i + are beyond the reach. 

Finally, we would like to emphasize that two quasi-degenerate heavy neutrinos of 
(9(100) MeV — O(10) GeV can account for not only the neutrino masses in oscillation 
experiments but also the baryon asymmetry of the universe [U El [6] . The heavy neutrinos 
studied in this work are thus quite interesting targets to search for. Furthermore, Ref. [7] 
reports that the sterile neutrinos with mass ~ 200 MeV could facilitate the energy trans- 
port from the supernova core to the schock front, prompting a successful explosion. The 
needed mixing is either |0 T | 2 > 10~ 8 or 10~ 7 — 10 -8 for |0 M | 2 . Interestingly, T2K can 
probe latter case via the N —y n~e + v mode (see the right pannel of Fig. [H]). In addition, 
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heavy neutrinos with masses smaller than (9(100) MeV may give a significant effect on 
the neutrinoless double beta decays [39]. According to a rough estimation in Table HJ 
MiniBooNE and SciBooNE have comparable abilities to T2K so that these experiments 
equally have the chance to probe these interesting possibilities. Serious analyses by these 
collaborations may lead to the discovery of the heavy neutrinos and revolutionize neutrino 
physics. 
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A Details on the flux calculation 

In this appendix, we present technical details on the flux calculation discussed in Section [3J 
Let us start from the source term Eq. PJ, 



As discussed in Section^ the neutrino flux 0„ (E„) of Eq. (J5]) is obtained by integrating 
this over 9, <p, I. The computation is greatly simplified if the kaon momentum is parallel 
to the beam axis. The detector geometry for this case is presented in Fig. [TTJ Under this 
simplification, the neutrino flux follows 

= ^ J^dl j l dcosQ S V {E V , 9, 1) P'(9, 9 ), (14) 

where A0 = r/i? sin #o and If = 96m. The function P'(9,9o) projects out the general 
polar angles into the ones with which pass through the near detector. It is given by 




(13) 



P'(6, 9 ) = H(9 - 9~)H{9 + - 9) 



(15) 



where H(x) is the Heaviside step function and 
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Figure 11: Geometry of the near detector ND280 and the kaon momentum. 



9(1) 



arcsm 



R sin #0 



R(l) = ^R 2 -2Rl cos9 + l 2 , 



(17) 



R(l) 

with R = 280 m and r = 3.5 m. 

The K + spectrum da/dpx in the proton collision with a graphite target is measured 
by NA61/SHINE Collaboration [30]. Fig. [TJ shows the best-fit values of the data [3D] and 
an example of the fit. The spectrum is the average over the K + polar angles (relative to 
the beam axis) 20 — 140 mrad where the probability that daughter neutrinos pass through 
ND280 is high [291 ED]- We assume that the shape of the kaon spectrum 4>k{pk) is not 
far from this measured spectrum and use the following fitting; 



k) 



fl 



( da'\ 
\dp K J 



da' 



?l<Ph 



a L p K 



I'L 



'>h = a H (p K +p ) bH , 
-(6 2 +i) 



;is) 



dp K 4>L + 4>H 

where a L = 1.285 mb GeV" (6l+1) , b L = 0.6289, a H = 793.5 mb GeV" l ° 2+i; , b H = 3.230, 
with which da'/dp K = da/dpx at p = 0. A positive value of p shifts the peak of da/dp K 
to lower energies. 

The differential decay width in Eq. ( 1T3"|) is given by 

d 3 T M 2 E u 



dE u dcos6d<f) 2E k 8tt 2 
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f(PK) 



where 



M 2 = 2G 2 F f 2 K m%\V u , 
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f(p K ) = m 2 K -m 2 - 2E v Jp 2 K + m 2 K + 2p K E v cos 9. 
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Figure 12: The K + spectrum da/dpx in proton-carbon collisions at 31 GeV proton energy. 
The dots are the best-fit values of the data 1301 and the solid line is the fit. 



An explicite expression of Eq. f JT4|) is written as 
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for < E v < (m 2 K — m 2 )/2m K and 
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for E v > (m 2 K — m 2 )/2mK- Here 
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For the heavy neutrino production K + — > fi + N, the differential decay width is given 



by 



where 



d 3 T 1 



dEjydcos 



|e M | 2 M x 2 PN 

2E K 8vr 2 
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2E N J p 2 K + m 2 K + 2p K p N cos 9. 



(2? 



(29) 



(30) 



The formula for K + — > e + N is obtained by making the replacement [i — > e in the above. 
By replacing the differential decay width in Eq. (Tl3l) with Eq. (1281) . it is straightforward 
to obtain the formulas for (J)n(En) similar to Eq. f l2"2"|) and Eq. 
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